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Q . Abstract 

Y^ • In this paper we introduce a new quantum computation model, the linear quantum cellular 

automaton. Well-formedness is an essential property for any quantum computing device since it 



enables us to define the probability of a configuration in an observation as the squared magnitude 
of its amplitude. We give an efficient algorithm which decides if a linear quantum cellular 

—^ ' automaton is well-formed. The complexity of the algorithm is 0(n-) in the algebraic model of 

^si , computation if the input automaton has continuous neighborhood. 
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0\ ■ 1 Introduction 

C^ ■ In order to analyze the complexity of algorithms, computer scientists usually choose some com- 

putational model, implement the algorithm on it and count the number of steps as a function 
of the size of the input. Different models, such as Turing machines (TM), random access ma- 
chines, circuits, or cellular automata can be used. They arc all universal in the sense that they 
can simulate each other with only a polynomial overhead. However, these models are based on 

C^ I classical physics, whereas physicists believe that the universe is better described by quantum 

mechanics. 

Feynman [^ |lj] and Bcnioff |Q, ^ were the first who pointed out that quantum physical 
systems are apparently difficult to simulate on classical computers, suggesting that there may 
be a gap between computational models based on classical physics and models based on quan- 
tum mechanics. Deutsch |lO[ introduced the first formal model of quantum computation, the 
quantum Turing machine (QTM). He also described a universal simulator for QTMs with an 
exponential overhead. More recently, Bernstein and Vazirani constructed a universal QTM with 
only a polynomial simulation overhead H . 

The power of QTMs was compared to that of classical probabilistic TMs in a sequence of 
papers jl^, |ll| |[ . The most striking evidence that QTMs can indeed be more powerful than 
probabilistic TMs was obtained by Shor|2^, who built his work on an earlier result of Simon 
PH . Shor has shown that the problems of computing the discrete logarithm and factoring can be 
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efficiently solved on a QTM, whereas no polynomial time algorithm is known for these problems 
on a probabilistic TM. 

Other quantum computational models were also studied. Yao p6[ has defined the quantum 
version of the Boolean circuit model, and has shown that QTMs working in polynomial time 
can be simulated by polynomial size quantum circuits. Also, physicists were interested in quan- 
tum cellular automata: Biafore M considered the problem of synchronization, Margolus |2G| 
described space-periodic quantum cellular automata and Lloyd llq, M] discussed the possibility 
to realize a special type of quantum linear cellular automaton (LQCA) . However these models 
are somehow different from the model of LQCA we consider in this article, and the physical 
realizability of our model has not yet been studied. 

Well-formedness is an essential notion in quantum computation. A quantum computational 
device is at any moment of its computation in a superposition of configurations, where each con- 
figuration has an associated complex amplitude. If the device is observed in some superposition 
of configurations then a configuration in the superposition will be chosen at random. The prob- 
ability a configuration will be chosen with is equal to the squared magnitude of its amplitude. 
Therefore it is essential that superpositions of unit norm be transformed into superpositions of 
unit norm, or equivalently, that the time evolution operator of the device preserve the norm. 
This property is called the well-formedness. In the case of a QTM, Bernstein and Vazirani gave 
easily checkable local constraints on the finite local transition function of the machine which 
were equivalent to its well-formedness. The existence of such relatively simple, local criteria is 
due to the local nature of the evolution of a TM: during a transition step only a fixed number 
of elements can be changed in a configuration. 

In this paper we will define formally linear quantum cellular automata and will give an 
efficient algorithm which decides if an LQCA is well-formed. Our algorithm is of complexity 
0{n^) if the input LQCA has continuous neighborhood (most papers in the literature in the 
classical context deal only with such automata). The problem of well-formedness in the case of 
an LQCA is much harder than in the case of a QTM. One cannot hope for local conditions on 
the local transition function as in the case of a QTM, since the transitions of a linear cellular 
automaton are global: a priori no constant bound can be given on the number of cells which 
are changing states in a step. It turns out that well-formedness is related to the reversibility of 
linear classical cellular automata. Thus our work is closely related to the decision procedure for 
reversible linear cellular automata of Sutner ES] . 

In fact, quantum mechanics imposes an even stronger constraint on any quantum compu- 
tational device: its time evolution operator has to be unitary. For QTMs ||j], space-periodic 
LQCAs [^ and partitioned LQCAs |25| well-formedness implies unitarity, but not for the model 
of LQCAs we consider here. Building on the present algorithm we gave in a subsequent paper 
III4 an efficient procedure which decides if the evolution operator of a LQCA is unitary. 



Watrous |25 has considered a subclass of LQCAs, partitioned linear quantum cellular au- 
tomata. He has shown that a QTM can be simulated by a machine from that class with constant 
slowdown, and conversely, a partitioned LQCA can be simulated by a QTM with linear slow- 
down. The efficient simulation of a general LQCA by a QTM is left open in his paper. As 
it is shown by Watrous, the problem of well-formedness in the case of a partitioned LQCA is 
easy. The local transition function of a partitioned LQCA can be described by a finite dimen- 
sional complex square matrix, and the automaton is well-formed if and only if this finite matrix 
preserves the norm. No analogous result is known in the case of a general LQCA. 

Our paper is organized as follows. In section 2 we first define linear cellular automata and 
give the basic notions of quantum computation in a finite space. Then we describe quantum 
linear cellular automata, define the notion of well-formedness, and prove that the inner product 
of two successor superpositions of configurations can be reduced to the inner product of two 
finite tensors. In section 3 first we give an example which shows that the trivial sufficient 
condition on the finite local transition function is not necessary for well-formedness. Then we 
describe the decision procedure for well-formed quantum linear cellular automata, prove its 
correctness, and analyze its complexity. The procedure consists of two separate algorithms, one 



which checks the unit norms, and another which checks the orthogonahty of the column vectors 
of the infinite dimensional time evolution matrix of the automaton. In section 4 wc describe 
a few open problems and finally in the appendix we give a shorter proof of one of the main 
theorems of Watrous' paper. 

2 The computation model 

2.1 Linear cellular automata 

A linear cellular automaton (LCA) is a 4-tuple A — (S, g, N ^ S). The cells of the automaton are 
organized in a line and are indexed by Z. E is a finite non-empty set of (cell-) states. At every 
step of the computation, each cell is in a particular state. The neighborhood N = (ai, . . . , a^) 
is a strictly increasing sequence of signed integers for some r > I, giving the addresses of 
the neighbors relative to each cell. This means that the neighbors of cell i are indexed by 
i + ai, . . . ,i + Or- We call r = \N\ the size of the neighborhood. Cells are simultaneously 
changing their states at each time step according to the states of their neighbors. This is 
described by the local transition function 5 : S'^' -^ S. If at a given step the neighbors of 
a cell are respectively in states xi,. . . ,Xr then at the next step the state of the cell will be 
6{xi, . . . ,Xr). The state g G S of ^ is the distinguished quiescent state, which satisfies by 
definition 5{q, . . . ,q) ~ q. 

The set of configurations is by definition S , where for every configuration c, and for every 
integer i, the state of the cell indexed by « is q. The support of a configuration c is supp{c) — 
{i G Ti : Ci y^ q}. A configuration c will be called finite if it has a finite support. We are dealing 
only with LCA's which work on finite configurations. Therefore from now on by configuration 
we will mean finite configuration. The set of configurations will be denoted Ca. 

The local transition function induces a global transition function, A : Ca ^ Ca, mapping 
a configuration to its successor. For every configuration c, and for every integer i, we have by 
definition 

[A{c)]{t)^S{c,+N), 

where 5(ci+Ar) is a short notation for ^{ci-^-al , ■ ■ ■ , Q+a^)- 

Configurations will often be represented by finite functions. We call an interval a finite 
subset of consecutive integers [j, k] ^ {j,j + I, . . . ,k} oi Z for any j and k (if j > k this defines 
the empty interval 0) . For our purposes it will be convenient to deal with representations whose 
domains are intervals. Therefore for a configuration c, and for an interval /, let c/ be the 
restriction of c to /. Also, let idom{c), the interval domain of c, be the smallest interval which 
contains supp{c). For an interval / = [j,k] with j < k, we define ext{I), the extension of / 
(with respect to the neighborhood N) as the interval [j — Or, k — ai]. The extension of is 0. 
If / == idom{c) then the support of its successor A(c) is contained in ext{I). Clearly, for every 
configuration c and intervals / and /', if idom{c) C / and ext{idom{c)) C /' then c/ and A(c)// 
specify respectively c and A(c). 

We will call an LCA simple if the elements of its neighborhood form an interval, that is 
a^ — ai = r — 1. In the literature LCA's are often by definition simple. 

A LCA is trivial if its neighborhood consist of a single cell. We can suppose without loss of 
generality that this single neighbor is the cell itself, that is A^ = (0). 

2.2 Basic notions of quantum computation 

Let _E be a finite set and let us consider the complex vector space C^ with the usual inner 
product which is defined for vectors u,v G C^ by 



The vectors in C^ will be called superpositions over E, and for a superposition u and an element 
e € E, we will say that u{e) is the amplitude of e in that superposition. The norm |[m|[ of a 
superposition u defined by this inner product is 



i«ii- JEk^)|' = v^ 



y eeE 

Two superpositions u and v are orthogonal, in notation m _L u, if {u,v) = 0. A superposition is 
valid if it has unit norm. If a valid superposition u over the set E is observed then one of the 
element of E will be chosen randomly and will be returned as the result of this observation. The 
probability that the element e is returned is |M(e)p. After the observation the superposition u 
is changed into the trivial superposition in which e has amplitude 1 and all the other elements 
0. 

Let / be an interval, and for each i g /, let Ui be a superposition over E. The tensor product 
®i^iUi is a superposition over E^ , that is an element of the complex vector space C^ , where 
by definition, for all x G E^ , 
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(x) = Y\_Ui{Xi). 

iei 

For our purposes the useful property of this operator is that the inner product of two tensors is 
the product of the respective inner products. Indeed, since I is finite, we have 
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=n^"-^*)- (1) 



2.3 Linear quantum cellular automata 

A linear quantum cellular automaton differs from a classical one in the sense that the automaton 
evolves on a superposition of configurations. The local transition function 5 maps the state vector 
of a neighborhood into a superposition of new states, giving the amplitude with which a cell 
moves into a specific state given the state of its neighbors. 

A linear quantum cellular automaton (LQCA) is a 4-tuple A — (E, g, N, 5), where the states 
set S and the neighborhood N are as before. It is called simple if the integers in TV form an 
interval. The local transition function is 5 : Tj'^' -^ C^ such that for every (xi, . . . , Xj.) G S*", 
we have ||5(a;i, . . . ,Xr)\\ > 0. The distinguished quiescent state g G S satisfies for all a; G S 
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Cells are again simultaneously changing their states at time steps but the outcome of the changes 
is not unique. If the neighbors of a cell are respectively in states xi, . . . ,Xr then at the next 
step, the cell will be in a superposition of states, where for every y G S, the state of the cell will 
be p with amplitude [(5(xi, . . . ,Xr)]{y). 

The local transition function induces a global one, which maps a superposition of config- 
urations into its successor superposition. We call it the linear time evolution operator Ua '■ 
Ca X Ca —> C For every c,d € Ca, the automaton enters d from c in one step with amplitude 

UA{d,c) = ll[5{c,+N)Kd,). 

iei 

This infinite product is well-defined since we deal with finite configurations, so for all but a finite 
number of integers i, q+at == (f and di ~ q. Therefore in the product only a finite number of 
terms can be different from 1. Moreover if there is an i such that c^+jv = Q^ and di ^ q then 
C/4((i, c) = 0. Thus in order to have non-zero transition amplitude it is necessary that idom{d) 
be contained in ext{idom{c)) . 



Let / be any interval which contains ext{idom{c)) . Then by the previous observations and 
by definition of tensor product we have 



UA{d,c) = 
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I (5(ci-f,Ar) [dj) if idom{d) C /, 

otherwise. 

Clearly, superpositions of configurations form the Hilbert space defined by 



HCa) ^{ueC^ 



y u{c) ■ u{c) < CXD > 
ceCA ) 



with the inner product defined for ui,U2 G C"-^ by 

(Ui, U2) = ^ -"1(0) • U2{c). 
c£Ca 

As usual, Ml and U2 are orthogonal (in notation ui 1. U2) if (mi,U2) — 0. 

As in the finite case, a superposition v of configurations is valid if ||w|| — \/ {v, v) = 1. 
Also, as in the finite case, if an LQCA is observed in a valid superposition of configurations v, 
the result of the observation will be the configuration c with probability |w(c)|^. Immediately 
after the observation whose outcome is c, the automaton will change its superposition into the 
classical one which gives amplitude 1 to c and to all the others. 

We want to have valid superpositions of configurations at each moment of the computation 
in order to associate the above probabilities to an observation. The initial configuration of 
the automaton is clearly valid. Therefore we say that the LQCA A is well-formed if its time 
evolution operator Ua preserves the norm. 

It is not hard to see that Ua preserves the norm if and only if its column vectors are orthonor- 
mal, that is they have unit norms and they are pairwise orthogonal. We will denote the column 
vector of index c by Ua{', c). In the next chapter we will give an algorithm which decides if the 
column vectors of Ua are orthonormal. An important technical tool in the correctness of the 
algorithm will be the generalization of equality ffl) to successor superpositions of configurations 
in the infinite Hilbert space. This is stated in the following lemma. 

Lemma 1 Let c and c' be configurations and let I be an interval such that ext[idom{c)) U 
ext{idom{d)) C /. Then we have 

{Ua{;c),Ua{;c')) 



n('5(Q+^),<5(c^+^)). 
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Proof 
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(2) 
(3) 

(4) 

(5) 
(6) 



The equations are justified in the following manner: (0) by definition of the inner product, (0) 
by the choice of /, (||) by identification of dj with d', (^ by definition of the tensor product and 
(|) by equation (0). D 

We have the immediate corollary: 

Corollary 1 Let c be a configuration and let I be an interval such that 
ext{idom{c)) C /. Then we have 

\\Ua{;c)\\=Y[\\5{c,+^)\\. 

3 A decision procedure for well-formed LQCAs 

3.1 Trivial LQCAs 

It is easy to give sufficient and necessary conditions for the well-formedness of a trivial LQCA 
which are easily checkable on the local transition function. 

Lemma 2 Let A = (S,g, (0),(5) be a trivial LQCA. Then A is well-formed if and only if for 
every a;, y G S with x ^ y 

Six) ± 6iy), (7) 

and for every x £ T. 

\\Six)\\ = l. (8) 

Proof For every a; e S let c^ be the configuration which is x at cell and quiescent elsewhere. 
Then for every x,j/ G S we have {S{x),S{y)) = {Ua{-,c^),Ua{-,c'^))- Thus if A is well-formed 
conditions (0) and (g) hold. 

For the converse suppose that both conditions are satisfied. Then corollary [| implies that 
the columns of Ua have unit norm. Now we show that for any two distinct configurations c and 
c', the associated columns of the evolution operator are orthogonal. Since c and c' are different 
there exist a cell i, such that a ^ c^. Thus d{ci) ± S{c[) by condition (0) and Ua{-, c) ± Ua{-, c') 
by lemma |l|. □ 

For non-trivial LQCAs condition (|^) can never hold since when |A^| > 1 we can not have 
|E|I''^I independent vectors in a space of dimension \T,\. 

But condition (g) still implies that the column vectors have unit norm by corollary ^. The 
following example shows that this condition is not necessary. 

Let B — {{q,p}, q, (0, 1), 6) be an LQCA with the local transition function defined as follows. 
For X G {q,p}, we define the superposition \x) over {q,p} by 

k)(y) = ( I 'J ^:^' 

' ''•'^^ ]^ if X j^ y. 
Then S is defined as: 

Siq,q) = \q), 5{q,p) = i|g), 

S{P,q) = 2|p), 5{p,p) = \p). 

In every configuration the number of pairs qp is equal to the number of pairs pq, therefore for 
all configurations c, d we have 

UBid,c) = l^ I If ^-J 

Thus the time evolution matrix Ub is just the identity, and B is well- formed. However, S{q,p) 
and S{p, q) do not have unit norm. 

6 



Nevertheless we can always transform a well-formed LQCA A — {Q, q, N, S) into an LQCA 
^' = {QtQi^i^') such that Ua — Ua' and A' satisfies condition (||). We simply renormalize 
the local transition function for all w S S'^' by defining 5'{w) — 5{'w)/\\5{w)\\. Then for every 
configurations c, d and interval / containing ext(idom{c)) and idom{d) we have 



UA'id,c) 



)(5'(ci+7v) 



lei 



idi) 



= Y[[S'{c,+N)Kd,) 



[S{ci+N)]idt) 

¥{C^+N)\\ 



- n 

lei 

U^eI\\Sic^+N)\\ 

\\Ua{;c)\\ 

1 

= UA{d,c). 

The following lemma establishes a particular property of trivial LQCAs which is not true in 
general. 

Lemma 3 Let A = (£,(?, (0),(5) be a trivial LQCA. If A is well-formed then Ua is unitary. 

Proof Suppose A is well-formed. By the previous lemma S is described by a unitary matrix. 
Let S~^ be the local function described by the inverse of this matrix, that is for all x,y G E 
we have [S-'^{y)]{x) = [5{x)]{y). Let A' be the trivial LQCA (E,g, (0),(5-i). Clearly Ua'Ua = 
UaUa' = I, which concludes the proof. □ 



3.2 The algorithm 

Before giving the algorithm, let us discuss the size of the input, that is the size of an LQCA 
A = {T,,q, N,5). It is clearly dominated by the size of the description of 6. We will work 
in the algebraic computational model, where by definition complex numbers take unit space, 
arithmetic operations and comparisons take unit time. Then S can be given by a table of size 
jEl''"''"'^, when the neighborhood is of size |iV| = r. Therefore we define the size of the automaton 
n = |S|''+^, and we will do the complexity analysis of our algorithm as a function of n. 
Our main theorem is an immediate consequence of Theorems and H. 

Theorem 2 There exists an algorithm P which takes a simple LQCA as input, and decides if 
it is well-formed. The complexity of the algorithm is 0{n^). 

What can we say about the well formcdncss of an LQCA which is not necessarily simple? 
Let A = (S, q, N, S) be an LQCA of size n whose neighborhood is A^ = (ai, . . . , ar). We can 
transform A into a simple LQCA A' — (S, q, N', 5') such that A and A' have the same time 
evolution operator. This can be done by taking as neighborhood iV' — (ai, oi-t-l, ai + 2, . . . , a^), 
and making the local transition function 5' independent from the new neighbors in N' . Then 
we can run P on A' . 

The size of A' will depend also on another parameter, on the span s of ^ which is defined 



as s 



ai 



1. Since |A'''| = s, the size of A' will be n' 



^s+ll 



j(a+l)/(r+l)^ Let 



define the expansion factor e of A as e = {s -\- l)/(r + 1). Then the time taken by P will be 
0(n'^) = 0{n^'^). We have therefore the following corollary: 
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Corollary 3 There exists an algorithm which takes an LQCA with expansion factor e as input, 
and decides if it is well-formed. The complexity of the algorithm is 0{n^'^). 

3.3 Unit norms of column vectors 

In this chapter wc wiU give an algorithm which decides if the column vectors of the time evolution 
operator have unit norms. Let A ~ (E, q, N, 6) be a simple LQCA whose neighborhood is of 
size r. We define an edge weighted directed de Bruijn graph Ga = {V, E, w) with vertex set 
V — S''^^, edge set E — {{xz,zy) : x,y € T,,z G E'""^} and with weight function w : E -^ R 
defined by w^^xz, zy)) — \\5{xzy)\\. The unweighted version of this graph was defined by Sutner 
in Ig^. A path is a sequence p = (woj • ■ • j I'fe) of vertices such that for 0<z<fc — 1, we have 
(vi, Vi+i) € E. The weight w{p) of a the path p is 

Y[ w{{v,,v^+i)). 

0<i<fe-l 

We call the path (wq, . . . , v^) a cycle if Vq = Vk and A: > 0. If in addition, vq — g''^^ then it is 
called a q- cycle. Our algorithm is based on the following lemma. 

Lemma 4 The column vectors of Ua have unit weight if and only if the weight of all q-cycles 
in Ga is 1. 

Proof Let T denote the set of g-cycles of Ga. We define a mapping M : Ca ^ T. Let c be 
a configuration with interval domain / — [j, k]. Let t = k — j, and for i = 0,1, . . . ,k — j, let 
Xi = Cj+i. Then by definition 

M{c) ^ {q''^'^,q''^^Xo,q''^^XoXi, ... ,XqX^ . . .Xr-2, ■■■ ,xtq''^'^ ,([^^). 

We have then 

||C/a(-,c)|| = \\5{q'-^x^)\\-\\8(<f-''x^xr)\\-...-\\8{xtq''-^)\\ bycoroUaryg 

= ||%OI|-|l%'^-'^o)||-...-P(x*g'-i)M|%OII 

— w{M{c)). by definition 

Since the mapping M is clearly surjective the statement of the lemma follows. □ 

Verifying if all column vectors of Ua are of unit norm is now reduced to checking if all 
g-cycles in Ga are of unit weight. The algorithm we give now will just do that. 

Theorem 4 There exists an algorithm R which takes a simple LQCA A = (E, q, N, S) as input, 
and decides if the column vectors of the time evolution operator Ua have all unit norm. The 
complexity of the algorithm is 0{n^). 

Proof Algorithm R will construct the graph Ga of lemma Q and then determines if it has a 
g-cycle of weight different from 1. This will be done by two consecutive algorithms i?i and i?2, 
from which the first will check if there is a column of norm less than 1, and the second will check 
if there is a column of norm greater than 1. They are both modifications of the Bellman- Ford 
single source shortest paths algorithm y, |l^, see also Q] (BF for short), when q^~^ is taken 
for the source. They are based on the fact that BF detects negative cycles going through the 
source. (Actually for our purposes any shortest paths algorithm can be used which uses sum 
and min as arithmetic operations, and which detects negative cycles. Floyd's algorithm would 
be another example). 

Algorithm Ri replaces every sum operation in BF by a product operation, and initializes the 
shortest path estimate for the source to 1 (the shortest path estimates for the other vertices are 
initialized to oo as in BF), and then runs it on Ga. This way it computes the shortest paths 
when the weight of a path is defined as the product of the edge weights. To see this let Ga' 



be the same graph as Ga except the edge weights are replaced by their logarithm. Then the 
weight of a shortest path in Ga given by BF will be the logarithm of the shortest path in Ga 
given by Ri . For the same reason, negative cycles in Ga through the source will correspond to 
g-cycles in Ga with weight less than 1 which will therefore be detected by i?i . 

Algorithm i?2 replaces every niin operation in Ri by max and the default initial shortest 
path estimate cxd by 0, and then runs it on Ga- This way it computes the shortest paths when 
the weight of a path is defined as the product of the reciprocal of the edge weights. If we define 
Ga with negative logarithm edge weights then negative cycles in Ga will correspond to cycles 
in Ga with weight greater than 1 and will be detected by i?2. 

The complexity of BF is 0{\V\ ■ \E\). In the graph Ga we have \V\ — |I]|'"^^. Every 
vertex has |E| outgoing edges, therefore \E\ — |S|''. Thus the complexity of the algorithm R is 
0(|Ep'-i) = 0(n2). D 

In |]l6| H0yer gave a linear time algorithm to decide if the column vectors have all unit norm, 
improving the complexity of our result. 

3.4 Orthogonality of column vectors 

Now we will build an algorithm which decides if the column vectors of the time evolution matrix 
are orthogonal. Let again A — (E, q, N, S) be a simple LQCA whose neighborhood is of size r. 
We define the graph Ha = {V,E) with vertex set V = E''^^ x E*"^^ and edge set 

E={ {{xiZi,X2Z2),{ziVl,Z2y2)) ■ 

xi,X2,Vi,y2 e E,zi,Z2 e Y7^'^,5{xiziyi) / 5{x2Z2V2) }• 

For a path p= {{uq^vo),. ..,{uk,Vk)) oiHA,\etpi = (mq, ■ • ■ ,Ufc), andp2 = (wo, • ■ •,«&)• Clearly, 
Pi and p2 are paths in Ga- A cycle is called here a q-cycle if its first vertex is (9'^^, q^~^)- 

Lemma 5 The column vectors o/Ua are orthogonal if and only if pi — P2 for every q-cycle p 
in Ha- 

Proof Let L — {{c,c') £ Ca ^ Ca ■ Ua{-tC) / Ua[-,c')}, and let T denote the set of g-cycles. 
We will define a mapping M : L ^ T. For (c, c') G L, let / = [j, k] be an interval such that 
ext(idom{c)) U ext{idom{c')) C /. Let t ~ k ~ j, and for z = 0, 1, . . . , A: — j we define Xi = Cj+i, 
and yi = c'_|_j. Then by definition 

M{c, c') = {{q^-\q'-'), {q'-'xo, q'-'yo), ■■■ , {x^q'-^ y^q'^-'), {q''~\q'-')). 

Since Ua{,-,c) / Ua{-,c'), lemma || implies that M{c,c') is indeed a g-cycle in Ha- Also, it is 
clear that M is surjective. Finally c ^ c' if and only if M(c, c')i 7^ M(c, c')2 since both are 
equivalent to the existence of i g / such that Xi ^ yi- □ 

We can now affirm: 

Theorem 5 There exists an algorithm S which takes a simple LQCA A — (T.,q, N,d) as in- 
put, and decides if the column vectors of the time evolution operator Ua are orthogonal. The 
complexity of the algorithm is 0{n^). 

Proof The algorithm S constructs the graph Ha and computes the strongly connected compo- 
nent of the node {q^~^, 'f~^)- By lemma rathere exists two distinct configurations such that the 
corresponding column vectors in Ua are not orthogonal if and only if in this component there 
is a vertex (w, v) with u ^ v. This can be checked easily. 

Finding the strongly connected components in a graph can be done in time 0(|-B|) for 
example with Tarjan's algorithm |2J]. In Ha the size of number of vertices is \V\ — jEp^'""^^. 
Since every vertex has outdegree |Ep, the number of edges is \E\ = jEp''. Therefore the 
complexity of the algorithm S is 0(|Ep'') = 0{v?). □ 



4 Conclusion 

It would be interesting to generalize results concerning reversibility of a linear classical CA 
for the well-formedness of an LQCA. For example a necessary condition for reversibility is the 
notion of balancedness of the local transition function H] , which means that every state has the 
same number of preimages. How does balancedness generalizes to the quantum model? 

It remains open, as stated also by Watrous, whether a QTM can simulate an LQCA with 
reasonable slowdown. 

Partitioned linear quantum cellular automata 

This appendix treats a special kind of LQCA, the partitioned LQCA, which was the main 
topic of Watrous' paper [Wat95]. Our aim is to provide a new, shorter proof to one of his results, 
based on our approach. 

A partitioned linear quantum cellular automaton (PLQCA) is a LQCA A = (S,g, A^, (5), 
which satisfies the following restrictions: 

1. The state-set E is the Cartesian product Ei x • • ■ x E^ of some finite non-empty sets E^, 
ie {l,...,r}. 

2. The local transition function (5 : E^ ^ C^ is the composition of two functions, the classical 
part (5p : E'' -^ E and the quantum part (5q : E — > C^. For all Xij G Ej, i,j G {1, . . . ,r}, 
Sp is defined by 

(5p((a;i,l, . . . ,Xi^r), (2^2,11 • • ■ ,X2,r), ■ ■ ■ , {Xr,l, ■ ■ ■ , Xr,r)) = (a^l,! , 2:2,2 , • ■ ■ ,Xr,r)- 

The function dp defines a LCA Ap = (E, N, 6p) whose global transition function Ap is a 
permutation on configurations such that for all c G Ca and i £ Z, 

[Ap{c)]{i) = Sp{ci+N)- 
Moreover, the time evolution operator Ua of Ap is a unitary matrix since for all c, d G Ca, we 

TT /^ ^^_ / 1 if Ap(c)=d, 
UA,(d,c)-<^ otherwise. 

The local transition matrix Q is the complex valued matrix, indexed by E, defined for all 
states X, y G E by 

Q{y,x) = [SQ{x)]{y). 

In fact Q completely determines the local transition function S. 

The function 5q defines a trivial LQCA Aq = (E, q, (0), (5q), with the time evolution operator 
Uaq ■ Clearly, Ca and q are respectively the set of configurations and the quiescent state also 
of Ap and Aq. It turns out that unitarity of the local transition matrix is equivalent to the 
unitarity of the time evolution operator, as stated in the following theorem. 

Theorem 6 ([Wat 95, theorem 3.1 and corollary 3.1]) Let A be a PLQCA, Ua its time 
evolution operator and Q its local transition matrix. Then the following statements are equiva- 
lent. 

L Q is unitary. 

2. A is well-formed. 

3. Ua is unitary. 

Proof The local transition function of A is the composition of two separate local transition 
functions, thus its time evolution operator is also the composition of time evolution operators of 
the associated LQCAs, that is Ua = UaqUa^,- Since Uap is unitary we have that Ua preserves 
the norm (resp. is unitary) if and only if UAn preserves the norm (resp. is unitary). 

The theorem follows from lemmas || and y. n 
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